Josephson Current in a Quantum Dot in the Kondo Regime Connected to Two 

Superconductors 
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We apply a Gutzwiller-like variational technique to study Josephson conduction across a quantum 
dot with an odd number of electrons connected to two superconducting leads. We show that, for 
small values of the superconducting gap, Kondo correlations and superconductivity cooperate to 
enhance the Josephson current. As the superconducting gap increases, the current changes sign and 
the system becomes a 7r-junction. The 7r-junction behavior sets in much before antiferromagnetic 
correlations at the dot can be treated perturbatively. 
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Phase coherent transport in hybrid superconducting- 
semiconducting nanostructures is already extensively in- 
vestigated. In these devices, quite interesting and sur- 
prising features emerge, due to electron-electron interac- 
tion In this letter we study a Quantum Dot (QD) at 
Coulomb Blockade (CB) with an odd number of electrons 
N, connected to superconducting contacts 0. Although 
such a kind of system has not been realized yet, it is 
very likely that continuous improvement of nanofabrica- 
tion techniques will soon make it available. In partic- 
ular, we address the question whether the formation of 
a Kondo singlet at the dot competes or cooperates with 
BCS s-wave pairing in the leads. We show that, as long as 
A <C ksTK, where Tk is the Kondo temperature of the 
dot and A is the superconducting gap, superconductiv- 
ity does not compete with Kondo correlation, but rather 
favors it. The system behaves as a regular Josephson 
junction and the critical Josephson current Ij is even en- 
hanced by Kondo coupling. However, as A becomes com- 
parable with IzbTk, there is a crossover to a regime where 
superconductivity spoils Kondo correlations on the dot, 
characterized by tt- Josephson coupling across the dot @. 

Coulomb Blockade has been widely investigated in 
Quantum Dots with normal leads 0] . If the coupling with 
the contacts is weak and the charging energy is higher 
than the thermal activation energy, DC conduction is 
strongly dependent on the gate voltage V g . The DC con- 
ductance across the dot shows a sequence of peaks, cor- 
responding to resonances between the chemical potential 
of the leads \i and an energy level of the dot. As the dot 
is tuned in a "Coulomb Blockade" (CB) valley between 
two consecutive peaks, iV is fixed and the conductance is 
heavily suppressed. 

In Superconducting/Normal/Superconducting (SNS) 
structures, Cooper pair subgap tunneling current takes 
place via Andreev states localized in the normal region 
||. When the normal region is given by a QD, charge 
quantization at the dot is not spoiled as the contacts be- 



come superconducting || . The I — V curve of a dot at 
CB between two superconductors has been derived in [|| , 
showing a rather interesting interplay between multiple 
Andreev reflection and electronic interaction at the QD. 

Kondo conductance may be achieved in a QD with 
N odd within a CB valley by increasing the coupling 
between dot and leads, provided the temperature T < Tk 
0. At T < Tk, a strongly correlated state between 
dot and leads sets in and a resonance in the density of 
states of the QD opens up at fi. Correspondingly, DC 
conductance across the dot increases, until it eventually 
reaches the unitarity limit at T = ||||. Since charge 
dynamics is "frozen out" by CB, the QD in this regime 
is usually modelized as a magnetic impurity with total 
spin S — I . 

Magnetic impurities embedded in a bulk superconduc- 
tor are known to strongly influence the superconducting 
critical temperature. Adding impurity states at energies 
below the superconducting gap can even give raise to 
gapless superconductivity Jul. 

Tunneling across a magnetic impurity between super- 
conductors with on site Coulomb repulsion has been re- 
cently revisited If A > UbTk, the system is un- 
able to scale toward the strongly-coupled Kondo regime. 
Sub-gap Cooper pair is strongly suppressed by Coulomb 
repulsion, unless each co-tunneling step is accompanied 
by a spin flip at the impurity. It has been proposed that 
such a mechanism may reverse the sign of the Josephson 
current through the dot so that the Josephson energy is 
at a minimum for <p = tt, where ip is the phase difference 
between the order parameters of the two superconductors 
attached to the impurity (7r-junction) fl^| . 

We expect that fine tuning of the parameters of a S- 
QD-S device (for instance, one may drive the system to 
a level crossing by tuning an external applied magnetic 
field which, under appropriate conditions, does not af- 
fect Kondo conduction [jp||l4]|)]]) may allow for the in- 
vestigation of the full range of physical conditions, from 
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A <C ksTx when the system can flow toward the strongly 
coupled Kondo regime prior to the onset of superconduc- 
tivity in the leads, to A > ksTx, when perturbation 
theory holds. 

We study the T = case with A <C ksTx using a non- 
perturbative variational technique. This technique has 
already been applied to the case of a dot with normal 
contacts, and it has been shown to provide good qualita- 
tive results in the perturbative regime as well, where it 
reproduces the poor man's scaling equation [ fl5| |. 

To construct the trial state, we start from the state 
\tp,s), given by the product of the left (L) condensate 
times the right (R) condensate, with the two order pa- 
rameters having a phase difference ip, times the state of 
the dot: Is). 



\<p,s) = |BCS, L) x |BCS, R, 92) x \s) 



(1) 



Here s is the spin component along the quantization axis 
of the dot spin = h, at CB with odd N. 

The minimal model for the Kondo interaction between 
electrons localized on the QD and elctrons from the con- 
tacts is Hk — Jct(O) • Sd- The spin density operator of 
the delocalized electrons, (?(0), at the position of the dot 
x = along the vertical axis, is : 



3(0) 



L,q,cr 



-L,q' ,a' 



+ C R,q',a') 



We have used the fermion operators c 9J (cJ .) (j = L, R) 
in the plane wave representation to describe the contacts 
particles and V is the normalization volume of the leads. 
We take the symmetric case in which hybridization of the 
dot with the L and R contacts, T, is the same. Kondo 
coupling is antiferromagnetic (AF): J > 0. The inter- 
action term Hk can be obtained from a Schrieffer Wolff 
transformation as in the case of normal contacts. The 
superconductivity in the contacts does not affect this 
derivation, provided D 3> A, where D is the bandwidth 
of the itinerant electrons. 

The total Hamiltonian is H = Hs + Hk, where 
H s, defined in eq.(^|) below, is the Hamiltonian for the L 
and R superconducting contacts in the BC S approxima- 
tion. 

The correlated trial state is constructed by applying 
a Gutzwiller-like projector P g to \<p,s). As T — > 0, 
the system scales towards the strongly coupled -large J- 
regime. Correspondingly, P g gradually projects out the 
high-energy components of the trial state, so that even- 
tually only a localized spin singlet survives at the QD. 

The "projector" P g is defined as |l5|| : 



g(a(0)) 2 ~ 4S d ■ <r(0) 



g is a variational parameter, which ranges between g = 
and g = 4/3. When g = 0, we have P = l-4S d -a(0). P a 



fully projects out the high energy localized spin triplet at 
x = 0. As g varies from to 4/3, also the localized spin 
doublet state it increasingly projected out. Eventually, 
when g reaches the value g — 4/3 , only the localized spin 
singlet is left over. 

The trial state is defined as: 



\g,ip) = P g \<p,s) 



(2) 



The value of g{J) is determined by finding the minimum 
of the energy functional E[g, J, ip, A], defined as: 



E[g,J,tp,A] 



{9,<p\g,<p) 



(3) 



Eq. (||) can be calculated by first expressing the products 
PgHPg and {Pg) 2 in terms of the usual fermion quasi- 
particle operators aj tq (ctj ), q ) (j — L,R) which 

destroy (create) excitations on the BCS states of the L 
and R contact and then by normal-ordering the corre- 
sponding operator products. The operators are: 



*h<i ~ u q c j,iA v q 



/3j- q = u qCj ^ qd +v q e l ^c] q 



(4) 



u q and v q are the BCS coherence factors and <pR = (p 
while 4>l = 0. The Hamiltonian for the contacts Hs is 
conveniently expressed in terms of these operators as: 



Hs — E] 



BCS 



+ B *(4« a i.« +$,-«&.-«)■ ( 5 ) 



E 

q.j=L,R 



Here £"bcs is the total ground state energy of the con- 
densates and E q = 4- A 2 are the energies of the 

quasiparticle excitations with £ q = q 2 /2m — /i, (we put 
h = ks = 1 throughout the paper). 

The variation of the energy due to the AF cou- 
pling at the quantum dot in units of the bandwidth D, 
e[£,j,8,(p] = E[g, J, A, (p]/D, takes a simple form once 
expressed in terms of the parameters £ = 5 (l — §<?), 
j = 3J/4D and 5 = A/D: 
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5 2 (l + cos(p))/N 2 (0)\ 2 



(1 + e) + (e - 1)<5(1 + cos( V .))/(2JV(0)A) 



2(1 - g + e) [VTTP - S 2 (l + cos(y))/JV(0)A] 
(1 + e) + (e - 1)^(1 + cos(^))/(27V(0)A) 



(6) 



Here A is the BCS electron-electron interaction 
strength and N(0) is the normal phase density of states 
at the Fermi level for each spin polarization. The first 
term is the expectation value of Hk, while the second 
is the raise in the average value of the kinetic energy 
Hamiltonian, Hs, due to the formation of the singlet be- 
tween the QD and the environment. The value of £ m i n , at 
which the energy is at a minimum measures how much 
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higher-energy spin states are projected out: £ m ; n = 
corresponds to full projection of states different from a 
localized spin singlet at the impurity. 

At d = the strong coupling fixed point is j — > oo and 
?min(j) — > 0. In correspondence of the fixed point, there 
is a large decrease of the Kondo energy: 



K 2 — TT(e-p — • (7) 

In Fig.|l| we plot £ m in vs. <5 for various j. As 5 is 
increased at j large, £ m ; n moves from to finite values. 
As (5 increases, two regimes can be envisaged: 
o) A -C Tk ■ The Kondo correlated state is firmly 
established. 

The change of £ m i n to first order in 6 is given by: 
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&nin(*) = £mta(* = 0) 



(5(1 + cos tp) 
N(0jX 



J + i - VP + 1 



that shows that, for any value of 5, ^ m ; n (^) < £min(<5 = 0). 
We also see that, as tp changes from to tt, £ m i n increases. 
This clearly shows that the stable configuration for the 
system is at tp = 0, what has important consequences for 
the Josephson conduction, as we are going to show in the 
following. 

In Fig.^ we plot e/j vs. £ for various values of j at 5 = 
0.03 and ip = 0. As j increases, £ m ; n moves toward £ m ; n = 
0, which corresponds to the strongly coupled fixed point. 
From the inset in FigJ| we also see that, as 5 increases, 
the value of the energy at the minimum decreases. In this 
regime e*[j, 8,<p] < e*\j, 0, tp] that is, superconductivity 
favors strong Kondo correlations. 

Upon calculating e[£,j, S, ip] for £ = £min> we get the 
best estimate for e* [j, 5, tp] , the energy of the correlated 
state, to first order in 6: 



e* [j, 6,<p]=CIk[1 + B6(1 + cos tp)] 



(8) 



lmn / ]). From 



where B = [1 - (£ nin ) 2 ]/(A^(0)A[l + (£ 
Eq.(||) we clearly see that, for large values of j, 
e*[j, 5, ip] < e*[j, 6 = 0, tp], that is, superconductivity and 
Kondo effect cooperate. 

b) A < Tk •' Superconductivity competes with Kondo 
ordering . 

Expansion in powers of S is no longer meaningful. The 
energy gain due to superconductivity Gs ~ —^N(0)A 2 
becomes dominant with respect to the Kondo energy gain 
&k{8 = 0). Kondo correlations and superconductivity 
start to compete. The minimum moves from tp = to 
tp = 7r, as shown in FigJ3 (bottom panel). Indeed, in 
the bottom panel of Figfswe plot e*/j vs. 5 for tp = 
and ip = 7r and j = 8. At a critical value of <5, the two 
curves cross each other, corresponding to crossing of the 
minimum energy from ip — to ip = tt. 

There is a third regime, not covered by our approach: 

c) A > Tk Superconducting order is dominant. 



Here the Kondo interaction can be treated perturba- 
tively. This regime was studied in Ref. J3| using a NCA 
perturbative approach. The results in p] show a change 
of the sign of the Josephson current that is consistent 
with our results, as we discuss below. 

The Josephson current is given by: 

Ij = -2ede*\j,5,<p]/d<p 



In regime a), Ij is linear in jS: 

Ij = — 2eflK BS sin ip 



(9) 



Therefore, in this regime Kondo correlations strongly 
increase the Josephson current. The sign of the current is 
the conventional one for Josephson systems, that is, the 
device behaves as a diamagnetic junction. As 6 slightly 
increases, the dependence of Ij on ip develops non sinu- 
soidal terms, in a way similar to the result in Ref. PH| . 

As 5 increases, we enter a new regime in which super- 
conductivity competes with Kondo correlations. In this 
region the Josephson current decreases as 5 increases (see 
top panel of Fig.^). As 6 — 5 P , the energy minimum 
moves from tp = to ip — tt. Correspondingly, although 
the current is now linear in S, its sign is reversed, that is, 
the device now behaves as a tt- Josephson junction. These 
results appear to extend the analysis performed in Ref. 
H in the regime A > Tk to the complementary regime, 
A < Tk- 

We now provide the physical interpretation of our re- 
sults in terms of Andreev bound states in the normal 
region. In the nonintcracting limit, the superconducting 
gap at the Fermi level \x in an SNS sandwich generates 
two localized Andreev states in the normal region with 
energy within the gap. The one of the two states with 
positive (negative) energy is mostly particle (hole)-like. 
In this case, Josephson conduction across the normal re- 
gion involves both Andreev states. 

As interaction is turned on ( J > 0) and A <C Tk, pair 
tunneling takes place through the Kondo resonance at fi, 
whose width (~ Tk) is quite large. At the same time, 
dot's spin "disappears" , because of full compensation at 
the strongly coupled fixed point ||. Full screening does 
not allow for extra occupancy of the dot. Accordingly, 
we expect that the hole-like state will be just above n but 
the particle-like state will be pushed quite high in energy 
(the higher is the effective coupling, the higher will be 
the energy of the particle- like state) . 

When superconductivity starts to weaken Kondo corre- 
lations, ( regime b) ), the two states are expected to have 
comparable energy, until, at some point, they cross. Such 
a crossing, where the particle state energy moves below 
the hole state energy and the particle state becomes par- 
tially occupied, makes the minimum move from ip = to 
ip = tt. By increasing A even further, we reach the regime 
discussed in ||] , where the particle state has moved below 
fi and the minimum of the energy has moved to tp = tt 
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(see bottom panel of Fig.||). This is the perturbative 
regime, where the Josephson current takes an overall - 
sign and the system behaves as a 7r-junction ||. 
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FIG. 1. £ m in vs. 8 for various j. For small values of 5, fmin 
is lowered, showing that a small amount of superconductivity 
favors Kondo correlations. On the other hand, for large 8, 
£ abruptly increases, that is a signal of the loss of Kondo 
regime. 



signature of Kondo effect in the quantum dot between 
two superconductors. 
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3 5 P 
FIG. 3. Top panel: saturation current I/(jD) vs. 8. For 

small values of 8, I/(jD) increases linearly with 8. At large 

enough 8 it decreases and becomes at S — S p , where the 

minimum of the energy moves to ip = n. For 8 > S p I < 0, 

that is, the system has crossed over to a 7r-junction behavior. 

Bottom panel: e*/j vs. 8 for ip = 0, n. At S — S p the two 

curves cross and the minimum moves to <p = ir. 
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FIG. 2. e/j vs. £ for different j and for tp = 0. At the onset 
of Kondo regime (large j), the minimum of the curve moves 
toward £ = 0. Inset: e/j vs. £ for two different values of <5. 
A small superconducting gap appears to lower e/j, that is, a 
small 8 favors Kondo regime. 

In conclusion, we have generalized the variational ap- 
proach introduced in ref. Jl5[ to study Josephson con- 
duction through a quantum dot at Coulomb blockade 
connected to two superconducting leads. Within such an 
approach, we analyzed the region of parameters where 
the dot lays in the strongly coupled Kondo regime and 
showed that, in this regime, it behaves as a diamagnetic 
junction, with a corresponding "Kondo enhancement" in 
the current. Our results appear to be in agreement with 
earlier predictions on related systems fll"|| . Moreover, our 
formalism shows that the transition to 7r-junction regime 
takes place much before antiferromagnetic correlations 
at the dot can be treated perturbatively, that is, much 
before Kondo effect has been disrupted by superconduc- 
tivity. Finally, our results lead us to predict an enhance- 
ment in the saturation current in the strongly coupled 
regime, which should be experimentally detectable as a 
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